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X ' Abstract 



We present an analysis of the nucleon strange sea extracted from a global Parton Distribution 
Function fit including the neutrino and anti-neutrino dimuon data by the CCFR and NuTeV col- 
laborations, the inclusive charged lepton-nucleon Deep Inelastic Scattering and Drell-Yan data. 
The (anti-)neutrino induced dimuon analysis is constrained by the semi-leptonic charmed-hadron 
branching ratio = (8.8 + 0.5)%, determined from the inclusive charmed hadron measurements 
performed by the FNAL-E531 and CHORUS neutrino emulsion experiments. Our analysis yields 
a strange sea suppression factor k{Q 2 = 20 GeV 2 ) = 0.62 + 0.04, the most precise value available, 
an ^-distribution of total strange sea that is slightly softer than the non-strange sea, and an asym- 
metry between strange and anti-strange quark distributions consistent with zero (integrated over x 
it is equal to 0.0013 + 0.0009 at Q 2 = 20 GeV 2 ). 
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1 Introduction 



The strange quark (s) distribution in the nucleon is an important input for the QCD phenomenology 
since the contribution of the s-quarks to the hard cross sections is of the same order of magnitude 
as the non-strange quarks. The strange quark contribution is particularly important at small values 
of the parton momentum fractions x, where the quark distributions are dominated by the sea. In 
high-energy hadron collisions the region of x < 0.1 is crucial for the study of many processes and 
therefore an accurate determination of the strange sea is required for the interpretation of experi- 
mental data. For instance a small positive s — s asymmetry in the strange sea may help explain the 
anomaly in the weak mixing angle reported by the NuTeV experiment (1). Inclusive cross sections 
are not very sensitive to the strange sea, since in this case the complementary contributions from 
strange and non-strange distributions are strongly anti-correlated. The strange sea is best con- 
strained by the neutrino-nucleon deep-inelastic scattering (DIS) dimuon data. This process stems 
from the charged-current (CC) production of a charm quark, which semileptonically decays into a 
final state secondary muon. The charm quark production cross section involves terms proportional 
to both the strange and the non-strange quark distributions. However, the contributions from u- and 
J-quarks are suppressed by the small quark-mixing Cabibbo-Kobayashi-Maskawa (CKM) matrix 
elements. The most precise (anti-)neutrino dimuon data currently available are by the CCFR and 
NuTeV collaborations [|2]|3l|4]|. In this paper we describe a determination of the strange sea distri- 
butions from a global parton distribution function (PDF) fit to the hard scattering processes, such 
as the inclusive charged-leptons DIS and Drell-Yan data, with the inclusion of the important CCFR 
and NuTeV dimuon data. The analysis is performed in the next-to-next-to-leading-order (NNLO) 
QCD approximation for the PDF evolution and for the massless coefficient functions. The next-to- 
leading-order (NLO) QCD corrections to the CC heavy-quarks production cross section are taken 
into account. These corrections reduce theoretical uncertainties on the strange sea due to variations 
in the renormalization and factorization scale. 

The paper is organized as follows. Section 2 presents the theoretical framework for (anti- 
neutrino induced charm dimuon production. In Section 3 we discuss the result of our global fit 
and the dominant theoretical uncertainties in the extraction of (anti-)strange quark distributions. 
We also discuss the impact of the semileptonic charm quark branching ratio on the strange 
distributions and we present an updated value of this parameter. Comparisons of these results with 
the earlier determinations of the strange sea from the leading-order (LO) analysis of Ref. [3 ], the 
NLO analysis of Refs. [|2l|5l, and the NNLO analysis of Ref. (61 are presented. Section 4 outlines 
future improvements in the determination of the strange sea distributions. 



1 



2 Theoretical Framework 



The differential cross section for charm quark production in CC (anti-)neutrino DIS off nucleon or 
nuclear target can be written as: 
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where x, y, and Q 2 are common DIS variables, E is the (anti-)neutrino energy, Gf is the Fermi 
constant, M and Mw are the nucleon and W-boson masses, respectively, and F2,r,3 are the corre- 
sponding structure functions (SFs). The nuclear data are usually presented in terms of an isoscalar 
target nucleon, which is the average over proton and neutron targets. For an isoscalar nucleon, 
assuming the usual isospin relations between the proton and neutron quark distributions, we have 
in the LO QCD approximation: 
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where u, d, s are the light quark distributions in the proton, £ = x{\ + ml/Q 2 ) is the slow-rescaling 
variable appearing in the kinematics of 2 — > 2 parton scattering with one massive particle in the final 
state [7J, and m c is the charm quark mass. The values of the CKM matrix elements V cs = 0.97334 
and V c d = 0.2256 [[81 suggest that the strange quark contribution dominates the cross section of 

Eq. (OQ) at small x. The factorization scale jx is usually set to either Q or ^Q 2 + m 2 . The sensitivity 
to a particular choice of \x gives an idea about the impact of higher-order QCD corrections. In the 
NLO QCD approximation the structure functions of Eq. Q get an additional 0(a s ) contribution 
from the gluon-radiation and gluon-initiated processes 0. In Fig. \T\ we compare the structure 
functions for charm production calculated in the NLO and LO approximations. The magnitude of 
NLO corrections rises at small x, giving the largest effect in the case of xF-$. For realistic kinemat- 
ics, the NLO corrections to Eq. © substantially cancel out in the difference between neutrino and 
anti-neutrino cross sections. In practice higher-order QCD corrections affect mainly the C-even 
combination s + s. We calculate the QCD-evolution of PDFs in the NNLO approximation IfTOl . 
However, a fully consistent NNLO calculation of the structure functions in Eq. (QQ) is currently 
not possible, since the NNLO coefficient functions for charm quark production are not available. 
The contribution to NNLO corrections from the soft-gluon re- summation has been calculated in 
Ref. [1 1J and is significant only at large values of x. Therefore, we do not include these corrections 
in our analysis. In general, the NNLO corrections are expected to be small compared to the un- 
certainties of experimental data, as one can infer from the typical magnitude of NLO corrections. 
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Figure 1: Comparison between LO (dashed) and NLO (solid) QCD approximations for charm 
quark production structure functions. The calculation is performed for neutrino interactions on 
isoscalar nucleons. 

We do not consider power corrections to the SFs of Eq. ©. The target mass corrections of 
Ref. |[T2l are marginal in the region of x < 0.3 covered by the CCFR and NuTeV dimuon data. The 
dynamical high-twist contributions to the charm production SFs are unknown. We estimate their 
effect by applying a simple rescaling for the quark charge to the phenomenological twist-4 terms 
extracted from the inclusive vN cross-sections [[T3l[T4| . Following this procedure we find that the 
impact of these corrections is negligible. 

Data from the CCFR and NuTeV experiments were collected on iron target. We apply nuclear 
corrections to Eq. © using the calculation of Ref. |[T5l [T6l . This calculation takes into account 
a number of different effects including the Fermi motion and binding, neutron excess, nuclear 
shadowing, nuclear pion excess and the off-shell correction to bound nucleon SFs. The model of 
Ref. lfT5Tl provides a good description of the nuclear EMC effect as measured in charged-lepton 
DIS over a wide range of nuclear targets, from deuterium to lead. In Ref. |[T6l this approach was 
extended to describe the (anti-)neutrino interactions with nuclei. The model predicts that nuclear 
corrections to the neutrino-nucleon structure functions are different from those for charged-lepton 
interactions. Furthermore, nuclear effects for the case of (anti-)neutrino scattering depend on the 
SFs type (F2 vs. xFt,) and on the specific C-parity and isospin states. Fig. [2] shows the nuclear 
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Figure 2: Ratio of differential cross sections for iron and isoscalar nucleon, ppe> in neutrino (left 
panel) and anti-neutrino (right panel) interactions [TT5l[T6l . The solid (dashed) curve corresponds 
to Q 2 = 3 (10) GeV 2 . The inelasticity, y, is fixed at 0.5. 

corrections for neutrino and anti-neutrino differential cross sections calculated for an iron target. 

Electroweak corrections including the one-loop terms are calculated in Ref. [fTTl . within the 
framework of the parton model, in a factorized form. In this approach the initial quark mass 
singularities of the QED diagrams are subtracted within the MS scheme and included into the 
PDFs, which absorb all electroweak corrections. It is interesting to note that the electroweak and 
nuclear corrections are similar in magnitude in certain kinematic regions. Since the dimuon data 
released by the NuTeV and CCFR collaborations have already been corrected for electroweak 
effects according to an earlier calculation of Ref. |fT8l , we do not apply such corrections in our fit. 

In the LO the dimuon cross section is related to the corresponding cross section for charmed- 
quark production as: 



where fh is the fraction of the charmed hadron h, D^(z) is the fragmentation function of the charm 
quark into a given charmed hadron h = D°,D + ,D^,A^ carrying a fraction z of the charm quark 
momentum, and Br(h — » /uX) is the corresponding inclusive branching ratio for the muon decays 
(note: the normalization 2/h = !)■ I n tne NLO the coefficient functions entering the SFs calcula- 
tion depend, in general, on z as well. The charm fragmentation function D c (z) defines the energy 
of the outgoing charmed hadron and, in turn, of the secondary muon produced in the semileptonic 
decays. Typically, a minimal energy E® is required for the muons identified experimentally, in 
order to suppress the background from light-meson semileptonic decays. Assuming a universal 
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Figure 3: Pulls of CCFR and NuTeV dimuon data with respect to our fit (left panel: neutrino, right 
panel: anti-neutrino). The solid lines represent a ±lcr band for the fitted model. The dashed dots 
illustrate the impact of an (anti-) strange sea enhancement on the (anti-)neutrino cross sections at 
small x. 



D c (z) for all charmed hadrons and integrating over z, Eq. (O reads: 

da-y^Eji > £°) _ B Jo-charm (4) 

dxdy ^ ^ dxdy 

where is the acceptance correction accounting for the cut > E®, and = ^fhBr(h — > fiX) 
is the effective semileptonic branching ratio. We use the values of rj^ evaluated by the NuTeV 
and CCFR collaborations lfT9ll . which are based on the NLO calculations of Ref. Il20ll and on the 
Collins-Spiller [|2TT| fragmentation function. The parameter s c , which defines the shape of D(z) in 
the Collins-Spiller model, is fixed at 0.6. This value corresponds to the best fit value obtained in 
the NuTeV analysis of Ref. |fT9t 

The charmed fractions fh depend on the incoming neutrino energy. This fact can be explained 
by the contributions from quasi-elastic A c and diffractive production. Furthermore, the values 
of fh are different for neutrino and anti-neutrino beams since in the second case no quasi-elastic 
A c production is present, but the relative rate of diffractive D s production is about a factor of two 
larger. These two contributions are significant mainly at low energies and they would not affect the 
value of at E v > 40 GeV. Measurements of fh and B^ in neutrino interactions were performed 
by the E531 E21I231 and CHORUS 1124112511 experiments using the emulsion detection technique. 
A value of B^ = 9.19 + 0.94% was obtained in Ref. (231 by combining the E531 data on f h in the 
energy range E > 30 GeV, which is relevant for the analysis of the NuTeV and CCFR dimuon data, 
with the charmed-hadron semileptonic branching ratios. The dominant source of uncertainty in 
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this determination of B^ is related to the uncertainties on the charmed fractions f^. 

A complementary determination of can be obtained from an analysis of dimuon data by 
performing a simultaneous fit of B^ with other parameters [2] El. In such an approach the absolute 
value of the dimuon cross section cannot directly constrain the (anti-)strange quark sea. This 
contribution is rather defined by the <2 2 -slope of the cross section, which is sensitive to the parton 
distributions through the QCD evolution equations. For anti-neutrinos the slope is driven mainly by 
the anti- strange sea, with a small contribution from gluons coming from the NLO corrections. In 
the neutrino case the non-strange quarks contribute as well. Once the (anti-) strange distributions 
are constrained by the <2 2 -slopes, the parameter B^ is determined by the absolute value of the 
dimuon (anti-)neutrino cross section. The value of B^ obtained from this global fit can then be 
compared with the direct measurements from the emulsion experiments in order to check the self- 
consistency between the <2 2 -slope and the absolute normalization of dimuon data. 

3 Results 

3.1 Constraints from CCFR and NuTeV Dimuon Data 

We determine the strange sea distributions from a global PDF fit to the CCFR and NuTeV dimuon 
data of Refs. [|2l SJ, combined with the inclusive charged-leptons DIS and the Drell-Yan cross 
sections used in the earlier fit of Ref. If26ll . The x-dependence of the strange and anti-strange quark 
distributions is parametrized independently using a model similar to that used for other PDFs: 

tkx,Q 2 Q )= A s x\l-x)" )s (5) 

at the starting value of the QCD evolution Qq = 9 GeV 2 . This functional form is flexible enough to 
describe the data. We do not observe any significant improvement in the quality of our fit by adding 
a polynomial factor to Eq. ©. The low-x exponents a s and a s are assumed to be the same as the 
one for the non-strange sea, since the existing dimuon data are not sensitive to them. The remaining 
parameters in Eq. © are extracted simultaneously with the non-strange PDF parameters, which 
essentially coincide with the ones obtained in Ref. [|26l . 

Our results for the strange sea parameters are given in Tabled! The quoted uncertainties include 
both statistical and systematic uncertainties in the data and take into account correlations in the 
latter where available. We obtain values of x 2 of 63 and 38 for the CCFR and NuTeV data sets, 
which both have 89 data points. It must be noted that, due to statistical correlations between 
data points, the effective number of degrees of freedom for the NuTeV data is about 40, which is 
consistent with our^ 2 value. 

The ratio of CCFR and NuTeV data with respect to the fit model is given in Fig. [31 Data from 
both experiments are consistent and are in agreement with our fit in the whole kinematic range. 
Although the CCFR anti-neutrino data is higher than the model at small x, this discrepancy is 
within the uncertainties. The dashed curves in Fig. [3] illustrate the effect of increasing the strange 
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Parameter 


Free B u 


Constrained B u 


A s 


0.086 + 0.007 


0.088+0.005 


a s 


-0.220 + 0.004 


-0.220 + 0.004 


b s 


7.7 + 1.0 


7.5+0.5 


A s 


0.083 + 0.008 


0.085 + 0.006 


a s 


-0.220 + 0.004 


-0.220 + 0.004 


b s 


8.0 + 0.4 


7.9 + 0.4 


m c (GeV) 


1.31+0.11 


1.32 + 0.11 


B„ (%) 


9.1 + 1.0 


8.80 + 0.45 



Table 1: Our results for the strange sea and charm production parameters. Central column: variant 
of the fit in which is extracted from the CCFR and NuTeV dimuon data only; right column: 
variant of the fit with B^ constrained by emulsion experiments. 



sea normalization parameters A s and A s by 0.03 @. One can see that both the normalization and 
the <2 2 -slope of the fitted model change with the strange sea normalization. If we model the 
energy dependence of B M by a linear function, the corresponding slope obtained from the fit is 
comparable to zero within uncertainties. We also do not observe any significant difference between 
the values of B^ obtained independently from the neutrino and anti-neutrino data sets: 9.4+ 1.1% 
and 8.9 + 2.2%, respectively. The neutrino-antineutrino and energy-averaged value of B^ = 9.1 + 
1.0 % obtained in our fit is in good agreement with the results of Ref. 11231 . We also extract the 
charm quark mass m c from the data. We obtain a value m c = 1 .3 1 +0. 1 1 GeV, which is in agreement 
with the world average MS value m c = 1.27+JJJ GeV 0. 

The strange sea suppression factor 

£ x\s(x,Q 2 ) + s(x,Q 2 )]dx 

<Q) = ~ = " — « ( 6 ) 

f Q x[u(x,Q 2 ) + d(x,Q 2 )]dx 

calculated with the PDFs obtained from our fit is given in Fig. HI The momenta carried by all sea 
quark flavors rise in the same way with Q 2 , due to the QCD evolution. Therefore, the suppression 
factor k also increases with Q 2 . We obtain k(20 GeV 2 ) = 0.59 + 0.08. The uncertainty in the strange 
sea normalization parameters is correlated with the one on B^. If we fix B^ at the central value 
obtained in our fit, we observe a reduction by a factor of 3 in the uncertainty on k. Our value of k 
is bigger than that obtained in the NLO QCD analysis of the CCFR dimuon data 0, k(20 GeV 2 ) = 
0.48+0 05. This difference occurs since the non-strange sea quark distributions used in Ref. (2) are 
larger than those of both our fit and other modern sets of PDFs (Fig. [5]). However, the strange sea 
from our fit is consistent with that of Ref. 0. The values of k calculated using the CTEQ6 11271 
and MSTW06 l|28l PDF sets agree with our determination within the uncertainties. The value 
of k preferred by the combined data on the vector meson electro-production in the analysis of 
Ref. ||29l is also consistent with our determination. The strange sea distribution obtained in our fit is 

4 We choose a shift corresponding to several standard deviations for illustration purpose. 
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Figure 4: Left panel: The +l<x band for the strange sea suppression factor k obtained in our fit 
(solid lines) compared to the determinations by the MSTW (dashes) and CTEQ (dashed dots) 
collaborations. The dotted lines represent the corresponding band after fixing the value of to 
the central value obtained in our fit, 9.1%. Right panel: The ±lcr band for the C-even combination 
of the strange sea distributions determined in our fit (solid lines) compared to the non-strange one 
scaled by k (dashes). 

somewhat softer than the non-strange one (see Fig.|4]). Due to the NLO corrections to the charmed- 
quark production coefficient functions the strange sea distributions are enhanced at small x (Fig.©. 
If we do not take into account such corrections, we obtain a smaller value of /c(20 GeV 2 ) = 0.55 ± 
0.13. This effect is consistent with the difference between the values of k obtained in the NLO 
fit of Ref . [O and in the LO fit of Ref . [3] . The variation of the strange sea due to a change of 



the QCD scale // from yJQ 2 + ml to Q is smaller than the one due to the NLO correction to the 
charmed-quark production coefficient functions. This result indicates our fit is stable with respect 
to the higher-order QCD corrections. 

In a variant of the fit with only the NuTeV data the strange sea is somewhat enhanced with 
respect to the combined CCFR and NuTeV fit (see Fig. [7]). The value = 7.2 ± 1.7 % obtained 
in this case is correspondingly smaller than those from both the combined fit and the analysis of 
Ref. [23 J. Although the discrepancy is at the level of lcr, it might indicate a certain inconsistency 
between the <2 2 -slope and the absolute normalization of the NuTeV data. In a variant of the fit with 
only CCFR data we get B^ = 9.7 +1.1 %, which is more consistent with the results from emulsion 
experiments. The strange sea determined from the CCFR data is somewhat smaller than the one 
from the combined fit. The strange sea charge asymmetry preferred by the NuTeV data is positive 
at all x values and is consistent with the analysis by the NuTeV collaboration @J. However, the 
CCFR data prefer negative charge asymmetry, so that the combined CCFR and NuTeV value is 
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Figure 5: The non-strange sea distribution obtained in the NLO QCD fit of Ref. @ by the CCFR 
collaboration (solid line) in comparison with the ones from MSTW06 (dashes) and CTEQ6 (dots) 
PDFs. The points give the corresponding strange-sea distribution extracted by CCFR 10. 

consistent with zero at the initial scale Qq = 9 GeV 2 (Fig. [7]). Once we impose the constraint 
s(x) = ~s(x), we observe an increase of x 1 limited to about one unit. The variant of fit with the 
constraint J Q [s(x) - s(x)]dx = imposed also does not yield statistically significant increase in the 
value of^ 2 . 

The strange sea asymmetry rises with Q 2 OOll because of the NNLO corrections. However, 
even taking into account such an effect, it remains consistent with zero within uncertainties in a 
wide range of Q 2 . In particular, at the reference scale Q 2 = 20 GeV 2 we obtain S~ = f Q x[s(x) - 
~s(x)]dx = 0.0010(13). The value of S ~ is sensitive to B^: if we fix the uncertainty on 5~ 
is reduced by about a factor of 2. The choice of the QCD scale fj. and the details of the high- 
order QCD corrections for the non- strange quark contributions to the charm SFs also affect S~ 

(Fig.©. Changing the QCD scale /u from ^Q 2 + m 2 to Q leads to a significant enhancement of the 
strange-anti-strange asymmetry at x ~ 0.15. The NNLO corrections to the QCD evolution and to 
the massless coefficient functions change the <2 2 -slope of the neutrino-nucleon DIS cross section. 
As a result, the strange sea distributions extracted from the fit, which are sensitive to this slope, 
are modified and the value of the strange asymmetry decreases. The nuclear corrections, discussed 
in Section 2, further reduce the asymmetry at x ~ 0.1. Each of these factors change the value of 
the asymmetry within 0.5cr. A combination of the effects discussed can, in principle, explain the 
difference between our result and those of Refs.|[51 0, in which a positive s — s asymmetry at the 
level of l-2cr was reported. 
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Q 2 =9 GeV 2 




x x 

Figure 6: Sensitivity of the strange sea distribution to various corrections and settings of the fit. 
Left panel: The shifts in C-even s + s distribution due to the NLO QCD corrections to the charm 
quark production coefficient functions (dashed-dotted curve), the variation of the QCD scale /u from 

sJq 2 + m 2 to Q 2 (dots), the NNLO corrections to the QCD evolution and the massless coefficient 
functions (long dashes), and the nuclear corrections (short dashes). The solid lines give the ±l<x 
uncertainty band from our fit. Right panel: The same for the C-odd distribution s — s. 

3.2 Impact of E531 and CHORUS Emulsion Data 



As explained in Section 13.11 the uncertainty on the strange sea derived from the fit can be sup- 
pressed if an additional constraint on the effective semileptonic branching ratio is imposed. 
Such a constrain can come from a direct detection of the charmed hadrons in the emulsion exper- 
iments. The only existing measurement of the charmed fractions fh as a function of the neutrino 
energy comes from a re-analysis ||23l of the data from the E531 experiment Il22l |3"TTl . Assuming 
/i-e universality and the recent values (8]| of exclusive branching ratios for charmed hadrons we 
can determine at different neutrino energies. Our results for the E531 data listed in Tables [2] 
and E correspond to B^D ) = 6.53 + 0.17%, B^{D + ) = 16.13 ±0.38%, B M (D+) = 8.06 ±0.76% 
and 5^(A+) = 4.50 ± 1.70% (SI and take into account correlations among the measured charmed 
fractions [|23l . Table [3] clearly shows that B^ increases with energy, with more pronounced varia- 
tions below 40 GeV. As explained in Section [2] the large contributions from quasi-elastic A c and 
diffractive production at low energies explain such energy dependence. Potential differences 
between neutrinos and anti-neutrinos are also expected to affect mainly the region E v < 40 GeV. 
This behaviour is consistent with the results of our fit to CCFR and NuTeV dimuon data described 
in Section I3T1 

The CHORUS experiment also measured the production rates of charmed hadrons in nuclear 
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Figure 7: Left panel: ±lcr bands for the C-even strange distribution s + s as obtained from the 
combined CCFR and NuTeV data (solid lines), from the CCFR data only (dashes), and from the 
NuTeV data only (dashed-dotted curves). Right panel: The same for the C-odd distribution s-s. 



Measurement 


E v > 5 GeV 


E v > 30 GeV 


CHORUS direct El 


7.30 + 0.82 


8.50+1.08 


CHORUS charmed fractions J23 


9.11+0.93 




E531 charmed fractions jj23j 


7.86 + 0.49 


8.86 + 0.57 


Weighted average 


7.94 + 0.38 


8.78 + 0.50 



Table 2: Semileptonic branching ratio B^(%) from direct measurements in the E53 1 and CHORUS 
emulsion experiments. The last row corresponds to our weighted average. 



emulsions. Thanks to a charm statistics about 20 times higher than the one of the E53 1 experiment, 
it was possible to directly detect some of the charmed-hadrons muon decays 112411 . The value of 
B M measured in Ref. [24] is given in the first line of Tabled A second independent measurement 
of can be obtained by combining the inclusive charmed fractions measured in Ref. 11251 with 
the corresponding branching ratios (HI, as explained above. The result is somewhat larger than the 
direct measurement as can be seen from the second line of Table [2] 

It is worth noting that all the determinations of B^ from emulsion experiments are sensitive 
to the value of the undetectable branching ratio D° — » all neutrals (0-prongs) 11241 [321 . which 
is decreasing the overall detection efficiency. The recent value for the fraction of 0-prong D° 
decays is 15 ±6% flU, which is intermediate between the ones assumed by the E531 and CHORUS 
analyses. 
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Figure 8: Left panel: The + l<x band for the ratio of the integral dimuon to the inclusive CC vFe 
cross sections as a function of the neutrino energy calculated using the results of our fit (solid 
curves). The charm production cross section ratio rescaled by the value = 8.8% is also given 
for comparison (dashes). Right panel: The same for dcr/ds, integrated over the neutrino energy 
spectrum of the NOMAD experiment [[37]| in the range of 6 + 300 GeV. A cut Q 2 > 1 GeV 2 is 
imposed in both cases. 



Energy (GeV) 


5 < E v < 20 


20 < E v < 40 


40 < E v < 80 


Ey >80 


B„ (%) 


6.33 + 1.05 


7.46 + 0.80 


8.68+0.85 


9.16 + 1.33 



Table 3: Semileptonic branching ratio for different neutrino energies obtained from the E531 
data ll23l and the recent values of inclusive leptonic branching ratios for D Q ,D + ,D + S ,K+ [[8]|. 



We can then proceed and average all emulsion measurements (see Table [2]). Uncertainties on 
such averaged values of B^ are smaller than the ones obtained from our fit to the CCFR and NuTeV 
dimuon data. The strange sea normalization is sensitive to variations of B^, so that the inclusion 
of the emulsion data on B^ to the fit reduces the uncertainties on the strange sea parameters. Since 
the energy dependence of B^ is more pronounced at small energies we use a single constraint 
Bp = 8.78 + 0.50% for E v > 30 GeV, as an additional data point in our global fit. Our independent 
extraction of B^ from the CCFR and NuTeV dimuon data, B fJ =9.\±\ .0%, is consistent with such 
measurement. Therefore, the central value of the strange sea parameters obtained in this extended 
fit are comparable with those obtained if B^ is unconstrained. However, the corresponding uncer- 
tainties are significantly reduced, as it can be seen from the second column of Table [TJ The value 
of the strange suppression factor becomes k(20 GeV 2 ) = 0.62 + 0.04, with an uncertainty twice 
smaller as compared to the variant of the fit with B^ unconstrained. With the constraint on B^ we 
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obtain a strange sea asymmetry S = 0.0013(9). This value is slightly larger than that obtained in 
the unconstrained fit, but still not significantly different from zero. 

4 Summary and Outlook 

In summary, we perform a global PDF fit using charged-lepton DIS data on proton and deuteron, 
fixed-target proton-proton and proton-deuteron Drell-Yan data, and (anti-)neutrino induced dimuon 
production data from CCFR and NuTeV experiments. We extract simultaneously the strange sea 
distributions and the effective semileptonic branching ratio B^ for charmed hadrons. The value 
of obtained by our global fit is consistent with the direct measurements from the E531 and 
CHORUS emulsion experiments. The constraint on B^ from emulsion data allows a reduction of 
the uncertainties on the strange sea parameters by about a factor of two. In particular, we obtain 
the absolute normalization of the strange sea with a precision of 6%, which is the most precise 
determination available. The x-shape of total strange sea is somewhat softer than the non-strange 
sea and the asymmetry between strange and anti-strange quark distributions is consistent with zero 
within uncertainties. 

An additional constraint on the strange sea distributions can be obtained from the inclusive 
(anti-)neutrino CC differential cross section dcr^/dxdy. At small values of x the scattering off 
strange sea quarks gives a significant contribution to the inclusive cross section. Available cross 
section data come from the CHORUS (33, NuTeV OH, and NOMAD fl35j experiments. The 
impact of the inclusive viV cross sections by CHORUS on the strange sea distributions was recently 
studied in Ref. IT361 in the context of a global PDF fit to the DIS data, resulting in a value of the 
asymmetry S~ = -0.001 +0.04. The inclusive CHORUS data were also included in an extended 
low-2 2 variant of our global PDF fit lfI51H51 . 

We expect a further improvement from the forthcoming measurements of the charmed fractions 
and the inclusive charm production cross section by CHORUS Il25ll . A global analysis of existing 
data from E531 and CHORUS emulsion experiments will allow a determination of B^ at a few 
percent level, improving the current dominant source of uncertainty on strange sea distributions. 

Finally, a sample of about 15k neutrino-induced charm dimuon events is expected from the 
ongoing NOMAD analysis [|38l . These data were collected on an iron target with an average beam 
energy of 24 GeV, and correspond to about three times the NuTeV dimuon statistics. Systematic 
uncertainties are kept well below statistical uncertainties through the measurement of the ratio 
of dimuon to inclusive CC cross sections, = cr^/crcc, as a function of different kinematic 
variables. Fig. [8]shows a prediction for the NOMAD experiment based on our current results. Pre- 
liminary studies indicate that the inclusion of the NOMAD dimuon data in a global PDF fit would 
substantially reduce the uncertainties in the determination of the strange sea distribution. Further- 
more, an accurate measurement of as a function of the partonic center-of-mass energy squared 
s = Q 2 (l/x-l) close to the charm production threshold would allow an improved determination 
of the charm quark mass m c . 
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